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Abstract 


Malaria continues to pose a major public health challenge, especially in developing countries, 219 million cases of malaria were 
estimated in 89 countries. In this paper, a mathematical model using non-linear differential equations is formulated to describe the 
impact of hygiene on Malaria transmission dynamics, the model is analyzed. The model is divided into seven compartments which 
includes five human compartments namely; Unhygienic susceptible human population S,, Hygienic Susceptible Human population 
Sn, Unhygienic infected human population /, , hygienic infected human population J, and the Recovered Human population R, and 
the mosquito population is subdivided into susceptible mosquitoes S, and infected mosquitoes /,. The positivity of the solution 
shows that there exists a domain where the model is biologically meaningful and mathematically well-posed. The Disease-Free 
Equilibrium (DFE) point of the model is obtained, we compute the Basic Reproduction Number using the next generation method 
and established the condition for Local stability of the disease-free equilibrium, and we thereafter obtained the global stability of the 
disease-free equilibrium by constructing the Lyapunov function of the model system. Also, sensitivity analysis of the model system 
was carried out to identify the influence of the parameters on the Basic Reproduction Number, the result shows that the natural 
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death rate of the mosquitoes is most sensitive to the basic reproduction number. 
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1. INTRODUCTION 


Malaria is one of the infectious diseases with an 





adverse effect on the human population. The 
malaria parasite that lives part of its life in humans 





and the remaining in vectors is transmitted between 
human host and mosquito vector by the infected 
female Anopheles mosquitoes through bites. In rare 
cases, people may be infected via contaminated 
blood, or a fetus may become infected by its mother 





during pregnancy or after delivery. Two of the five 
parasites species — Plasmodium Falciparum and 
Plasmodium vivax- pose the greatest public health 
challenges [1] — [3]. According to the World Health 
Organization [4], 219 million cases of malaria were 





estimated in 89 countries, the estimated deaths were 
435,000 with the 
disproportionately high share of the global malaria 


African region carrying a 


burden. Malaria is the third leading cause of death 
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most especially for children under five years, after 
pneumonia and diarrheal disease. 

One of the factors that increase the transmission 
dynamics of malaria in poor environmental 
conditions. These conditions include unsafe water 
supplies, poor personal hygiene, poor sanitary 
facilities, poor living standards, and unhygienic 
food. Poor personal hygiene may result in water- 
borne diseases [5]. According to Nkuo-Akenji et al. 
[6] poor environmental sanitation (hygiene) and 
housing conditions might be significant risk factors 
for malaria burden. Also, Enebeli et al., [7] in their 
study concluded that poor access to water, 


sanitation, and hygiene practices of caregivers 





directly relates to the prevalence of malaria among 
their children. Olaniyi et al. [3] developed a 
mathematical model of malaria dynamics with 
in the 


presence of protected travelers. Okuonghae and 


naturally acquired transient immunity 
Nwankwo [8] developed a non-autonomous model 
to assess the impact of different microclimate 
conditions on the transmission dynamics of malaria. 
Goni and Musa [9] proposed a mathematical model 
for the transmission dynamics of malaria by 
incorporating change via education as a control 
strategy. The human population follows the 
susceptible-protected-exposed-infectious-recovered 
(SPEIR) pattern and the mosquito population 
(SED) 


patterns. An analytical study was carried out to 


follows —susceptible-exposed-infectious 
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investigate the local stability of the system, the 
basic reproduction number was obtained using the 
next-generation matrix method. The result shows 





that the disease-free equilibrium of the system is 
locally asymptotically stable if Rə < 1. Bhuju et al. 


[10] mathematically studied the impact of 
temperature in malaria disease transmission 
dynamics. The SEIR model for the human 


population and LSEI compartment model for 
that 
temperature affects the transmission dynamics of 
malaria significantly. Okosun and Makinde [11] 


mosquito population. It was observed 


also modeled the impact of drug- resistance in 
malaria transmission. Many mathematical models 





have been developed to study malaria dynamics but 
none has been discussed to study the impact of 
hygiene on malaria transmission dynamics as 
proposed in this work. 

In this work, we propose a deterministic 
mathematical model for assessing the impact of 
hygiene on malaria transmission dynamics. We 
compute the Basic Reproduction Number and 
establish the local and global stability of the disease 
-free equilibrium and carry out the sensitivity 
analysis of the parameters. 


2. MATERIALS AND METHODS 


2.1. Model Formulation 

In this model, the total human population 
subdivided 
human population Sy, 


denoted by Ny is into Unhygienic 


susceptible Hygienic 
Susceptible Human population Sa, Unhygienic 
infected human population Iw, hygienic infected 
human population 7, and the Recovered Human 
population r,. The mosquito population denoted by 
N, is subdivided into susceptible mosquitoes S, and 
infected mosquitoes I,. So 


Ng = da + Sn t i. + In + R. (2.1) 


Ny = Sy + Ip. (2.2) 

Let Ay be the recruitment rate of the human 
population. A fraction (1—@)Ay enters unhygienic 
susceptible human class while the remaining 
fraction @Ay enters the hygienic susceptible human 
class. The unhygienic susceptible class 1s increased 
by the rate at which unhygienic human class lose 
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immunity after recovery given as w,, and reduced 
by the rate of progression to hygienic class qı, the 
force of infection for the unhygienic class 4, and 
natural human death rate my. The hygienic 
susceptible human compartment is increased by 
the tı, the rate at which hygienic human loss 
while the 


compartment is reduced by natural human death 


immunity after recovery at wp, 
rate # and the force of infection for the hygienic 
class (1—%4» The unhygienic infected human class 
I, is increased by 4, and reduced by natural 
human death rate ug, rate of progression from Iu to 
In given as t2, malaria induced death for unhygienic 
human class è, and recovery for unhygienic 
infected class J, is 
increased by (1—%)4 and tə then reduced by the 


recovery rate for a hygienic human class given 


human #@,. The hygienic 


as ĝa, malaria induced death for hygienic human 
class 6, and natural death rate ug. The Human 
recovery class R is increased by, and #,, then 
reduced by Hg, ©, and w,. The susceptible 
mosquito class S, is increased by the Mosquito 
A,, reduced by the 


mosquitoes death rate p, and force of infection for 


recruitment rate given as 


mosquito given as 4». The infected mosquito class 
I, is increased by A, and py. 


Given the above description and definitions of 
variables and parameters in Table 1 and 2, the 
following are the model equations: 





dSy | 
i (1-@)Ag — (Ti + Au + Hn)Su + OB, (2.3) 
dS 
= = Wp + WR + T1Su — (AU — DAn #p)Sp, (2.4) 
dl, 
ae AS — (t2 +6, +9, + Hay, (2.5) 
di 
= = (1 — An Sp -= Taly == (ôn + On + Ep), (2.6) 
dR | | 2 
R = Ouly + Onn — Wu + On + Hy )R, (2.7) 
dS, 
dt Ay — AySy — HySy, (2.8) 
dl, 
== AySy — Pyle (2.9) 
“where” 
_ Di Bvhlv _ baBohlv 
Au = a , Àn — Nu , b; > bz, 
L= ba Biv (tut eth) g > En, On >O. (2.10) 


Ny 
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Figure 1. Model flow diagram. 


2.2. Invariant Region 
The invariant region can be obtained by the 
t>0 following theorem. 


Theorem 1 
The solutions of the model are feasible for all 
if they enter the invariant region 


O = Oy XO). (2.11) 
Proof: 

Let 
0 = (Su Sn Iu, In R, Sy, Ip) ERI, (2.12) 


be any solution of the system with non-negative 
initial conditions. Hence, all feasible solution set of 
the human population of the malaria model enters 
the region 


Qy = (Su, Sn, L In, R) € RS: S, > 0,5, > 0, 
I, = 0, I, 2 0,R 2 0, Ng < a} (2.13) 
Similarly, the feasible solution set of the vector 
population enter the region 





O, = {(S,1) E RẸ: S, 2 0, ly 2 0, Ny <=}. (2.14) 


Therefore, the region Q is positively invariant 


i.e. the solution remains positive for all initial 
values. 

Thus, the model is biologically meaningful and 
mathematically well-posed in the domain 2. 


2.3. Disease Free Equilibrium (DFE) 

The DFE of the model equations can be found 
by setting the right hand of the model (2.3) — (2.9) 
to zero, 1.6. 


dS, dS, di, diy dR dS, di, 


dt dt dt dt dt dt dt ` 


which gives 





(1 — @)Ag — (7, +A, + Wg) S, + @,R = 0, (2.3) 
aAy + w,R + 7,5, —(G — DA, +Hy)Sh = 0, (2.4) 
Table 1. Variables. 
Symbols Description 
Sy Unhygienic Susceptible Human 
Sn Hygienic Susceptible Human 
Ia Unhygienic Infected Human 
I Hygienic Infected Human 
R Recovered Human 
Sy Susceptible Mosquitoes 
J, Infected Mosquitoes 
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Table 2. Model Parameters. 


Symbols Description 
Ag Recruitment rate of Human Population 
Ay Recruitment rate of mosquitoes 
Ti Progression from S, to Sn 
Tə Progression from l, to l, 
Ô Disease—induced death for the unhygienic human class 
On Disease—induced death for the hygienic human class 
b, Biting rate of mosguito for unhygienic human class 
b- Biting rate of mosquito for hygienic human class 
Bon Transmission probability of infection from mosquito to human 
Bry Transmission probability of infection from human to mosquitoes 
Ay The force of infection for unhygienic human class 
An The force of infection for hygienic human class 
Ap Force of infection for mosquitoes 
bz Biting rate of mosquitoes 
4 Rate of reduction of infection for hygienic class 
p Modification parameter 
a, Rate of recovery for unhygienic human class 
n Rate of recovery for hygienic human class 
W Rate at which recovered human become susceptible 
a Hygienic rate 
Ly Natural human death rate 
Lly Natural death rate of mosquitoes 
Ny Total human population 
AuSu — (12 + Ôu + Ou + Hey = 0, (2.5) At DFE, Iu = In = Iv = 0, 
(1-DAnSn + Talu — (őn + Oy + uy, = 0, (2.6) So we have, 
Oulu + OnIn — (Wu + wn + pp)R = 0, (2.7) (1 — @)Ay — (7, + Hg )Su = 0, 
Ay — AySy — HySy = 0, (2.8) aAg+t,S, — HySn = 0, 
AS, — ttyly = 0 (2.9) A, — HyS, = 0. 





IGD PANDAWA 4 


J. Multidiscip. Appl. Nat. Sci. 


After computing simultaneously, we have 


s = (1 —a@)Ay 
“ (+ ux)’ 


Ag(t, + apy) 
5, = — 
Hy (Ty + Hg) 


Therefore, the Disease Free Equilibrium (DFE) 
point of the model is given by 





(1-a)Ay Ag(t1+Guy) 
(Tithe) Matra) ' 





By (Sse ie St 
(2.15) 


2.4. Reproduction Number (Ro) 

The basic reproduction number (Ro) 1s defined 
as the number of secondary malaria infections 
produced by one infected individual in a completely 
susceptible The 
method [12] will be employed to compute Ro. F(x) 
is the rate of new infection appearance while V(x) 


community. next-generation 











is the rate of transfer of individuals into 
compartments. So we have 
(Ty +H) 
(1—€) bz vh (QU +14) 
— 0 0 RESA SZ (1G 
F (Ty +H) ( ) 
ba BhyAyiiq PBzbhrðvytH 0 
Ap ly Apgby 
k 0 0 
V = (~z ka 0 ) (2.17) 
0 0 m 
“where” 
ky = (2 + 6, + Ou + ug), kz = (On + On + Hp), 
kz = (© + py). (2.18) 
= 0 
Ky 
tr—l1 T3 1 
ia kiko ka 0 |, 
o 0 - 
By 
Thus, 
0 0 bi Bon (1-0) up 
(Tr t+HH) My 
e 0 0 -Dbz Bonlanatra) | 
(Ty +H) My 
bz BhyAyiiq TzPbz3BpyAyig  Pb3BhvAybH 0 
Apbyki AppykyKz Apbykz 
(2.19) 


0,0,0,°, 0). 


The basic reproduction number is the largest 
eigenvalue or spectral radius of FV’. Hence, 





R, - bs BvhEhvAvbn(bpxp(1-@)(k2+T2p)+b,k1p(aup+t)(1-%)) 
7 4 AUBky ko (ty +My) | 
(2.20) 
Theorem 2 
The equilibrium Eo (DFE) for 


model system is locally asymptotically stable if 





disease-free 
Ro < 1 and unstable otherwise. 


Proof: 
At DFE, the Jacobian matrix is given by 


—(t1 tp) 0 0 0 w 0 0 
Tı — Uy 0 0 Ww 0 
0 0 —k, 0 0 0 0 
ia 0 0 0 -k, 0 0 0 
0 0 0 0 -k 0 0 
0 0 0 0 0 uu, 0 
0 0 0 0 0 — [Ly 





Table 3. Indices of Sensitivity. 


Symbols Sensitivity Index 
Ay -1 
A, l 
Tj -0.00013 
Tə -0.000041 
On -0.000041 
On -0.29 
b, 0.00022 
b- l 


Bon l 


b3 l 
-0.087 

p l 

Ou -0.000015 
On -0.71 

a -0.00011 
HH l 

Hy -2 
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y — Atl 
—(tı +g) —A 0 0 0 w 0 0 
Ti —tn —A 0 0 w 0 0 
0 0 | 0 0 0 0 
= 0 0 0 = oe | 0 0 0 
0 0 0 0 or | 0 0 
0 0 0 0 0 —H, —A 0 
0 0 0 0 0 0 | 


The eigenvalues are: 


Ay = —(t + My), Az = — bg, Az = —ky, Ag = —ka, 
A, = —kz,Ag = Az = — Hp: 


It is observed that all the eigenvalues are 
negative, this implies Ro < 1 that at the DFE point 
is locally asymptotically stable, this means that 
Malaria infection can be eliminated from the 
population. 


Table 4. Parameter values of Model. 


Symbols Values 
Ag 100 
A, 1000 
Ti 0.25 
Tə 0.5 
Oa 0.13 
On 0.06 
b, 0.17 
bz 0.1 
Bon 0.03 
Bry 0.09 
bz 0.12 

0.08 
0.5 
Bb, 0.05 
n 0.15 
a) 0.7902 
a 0.46 
Hp 0.00004 
Ho 0.0000569 
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2.5. Global stability of the Disease Free 
Equilibrium (DFE) 
Theorem 3 

The DFE of the model system is globally 
asymptotically stable if Ro = 1. 


Proof: 
Consider the following Lyapunov function: 


V(t) = b3aBnoAphtg(k2 + T2P)Iu + bapApttaBnokiIn 


+ Agkykopyly, (2.21) 


Differentiating we have 


av 
a. bz Bro oli (ke + Tp) (A, Sy — kilu) + b3pAobgBnoki 
((1 — Q) Ap Sp + Talu —k2In) + Agky kop, AS, — Holy), 
(2.22) 


Source 
[13] 
[14] 

(Assumed) 

(Assumed) 

(Assumed) 

(Assumed) 

(Assumed) 

(Assumed) 
[2] 

[2] 
[15] 
(Assumed) 
(Assumed) 
(Assumed) 
(Assumed) 
[14] 
(Assumed) 
[13] 
[14] 
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At DFE, we have 


dt (Ti ttp) 


— Apkikotis ) Ty, 
dt | Ap ka ko 45 lTi ttp) 
— 1) L, 


dF 
— < (RẸ — Dip (2.23) 


From the equation above, < < o, if Ro = 1. (2.24) 


Hence, the DFE is globally asymptotically 
stable. 


2.6. Sensitivity Analysis 

In this section, sensitivity analysis is carried out 
to identify the parameters that have a great 
influence on the basic reproduction number (Rọ). 
The sensitivity index of Rp to a given parameter P is 
given by the relation 


ee (2.25) 
Table 3 shows the sensitivity indices of the basic 


The 
parameters with positive indices indicate that the 


reproduction number to the parameters. 


basic reproduction number increases as their values 
increase. While the parameters with negative 
sensitivity indices indicate an increase in these 
parameters will result in the decline of the basic 


reproduction number and vice-versa. 
3. RESULT AND DISCUSSIONS 


First, we show that this system of model is 
biologically meaningful and mathematically well- 
posed in the given domain 0. We compute the basic 
reproduction number Ry, of the model using the next 
-generation method. We established the existence of 





the disease-free equilibrium of the system and show 
the condition for the local stability of the disease- 
free equilibrium and global stability of the disease- 
free equilibrium using the Lyapunov function. The 





disease-free equilibrium is locally asymptotically 
stable 
if Ro = 1. Sensitivity analysis of the model equation 


if Ro < 1 and globally asymptotically stable 


was carried out as illustrated in Table 3. From the 
table, it shows that the natural death rate of 


sf, (Patton Provi Paru 0 +120)+b2k1p(aup+t,)(1-%) 


e ee 


mosquitoes ( H» ) is most sensitive to the Basic 
Reproduction Number. 


4. CONCLUSIONS 


In this work, the mathematical model to 


assess the impact of hygiene on malaria 
transmission dynamics was proposed and analyzed. 
The model is divided into the Human Population 
and Vector (Mosguito) Population, the Human 
population is further subdivided into the Susceptible 
Susceptible 


hygienic Human population, infected unhygienic 


Unhygienic human Population, 


human population, infected hygienic human 
population and Recovered human, while the vector 
population is subdivided into the susceptible vector 
and infected vector. We proved that the model 
meaningful and 
The 


equilibrium (DFE) was established and it was 


equation is biologically 





mathematically  well-posed. disease-free 


observed that DFE is locally asymptotically stable 
Ro<1if 
Ro = 1 using 


and globally asymptotically stable if 


Lyapunov function. Sensitivity 
analysis of the model parameters was carried out 
and it shows that the natural death rate of 


mosquitoes is most sensitive to the Basic 
reproduction Number. 

This implies that individuals must continue to 
engage in activities that promote both personal 
hygiene and environmental hygiene so as reduce the 
growth of mosquito hence curbing the spread of 
Non- 


Governmental Organizations (NGOs) must continue 


malaria also government and other 
to intensify campaigns on hygienic practices at 
individual and community levels. Future studies can 
be carried out on the model this includes: establish 
and prove the existence of the unique endemic 
equilibrium point, Stability analysis (local and 
global stability) of the endemic equilibrium point, 
and lastly solve the model equations using any 


analytic method available. 
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